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Aims. The properties of young stellar clusters (YSCs) in M33, identified from the center out to about twice the size of the bright 
star-forming disk, are investigated to determine possible spatial and time variations of the star formation process in this Local Group 
' blue galaxy. 

Methods. 915 MIR sources have been extracted from the Spitzer 24 fim image. Upon inspection of Ha and GALEX images and 
exclusion of evolved AGB stars, a sample of 648 objects is selected as candidate YSCs and their luminosity function is examined. The 
spectral energy distribution of each object, based on aperture photometry, is compared with Starburst99 models to derive age, mass 
and A v of individual clusters. In the analysis we allow for different values of the upper mass cutoff of the stellar initial mass function 
(IMF), the porosity of the ISM, and the dustiness of HII regions. We also examine the influence of different dust models and include 
corrections for incompleteness of the IMF. 

Results. We find discrete MIR sources as far as the extent of the warped HI disk, i.e. 16kpc from the galaxy center. Their surface 
density has a steep radial decline beyond 4.5 kpc, and flattens out beyond the optical radius at 8.5 kpc. We are able to identify YSCs out 
to 12 kpc. At large galactocentric radii, the paucity of luminous clusters and the relevance of hot dust emission become evident from 
the analysis of the bolometric and MIR luminosity functions. The YSC mass and size are correlated with a log-log slope of 2.09 + 0.01, 
similar to that measured for giant molecular clouds in M33 and the Milky Way, which represent the protocluster environment. Most 
of the YSCs in our sample have Ay ~0-l mag and ages between 3 and 10 Myr. In the inner regions of M33 the clusters span a wide 
^ ' range of mass (10 2 < M < 3 x 10 5 M Q ) and luminosity (10 38 < L bol < 3x 10 41 erg s -1 ), while at galactocentric radii larger than 

~ 4 kpc we find a deficiency of massive clusters. Beyond 7 kpc, where the Ha surface brightness drops significantly, the dominant 
YSC population has M < 10 3 M Q and a slightly older age (10 Myrs). This implies the occurrence of star formation events about 
lOMyr ago as far as 10-12kpc from the center of M33. The cluster L FUV vs. L Ha relation is non-linear for L FU v < 10 39 ergs~', in 
agreement with randomly sampled models of the IMF which, furthermore, shows no appreciable variation throughout the M33 disk. 
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^J. , 1 . Introduction physical scales which can be spatially resolved will never be as 

Q\ ■ small as in Galactic surveys, some global properties of the YSCs 

<3 Young stellar clusters (YSCs) are important tools for the study of can be investigated using their integrated spectral energy distri- 

i the star formation (SF) process in galaxies. Having just formed butiori5 fro m the UV to long enough wavelengths to overcome 

, from gravitationally bound molecular clouds and often still em- extinction dGrossi et al.l2010h . Tracers of SF such as IR, Ha, far- 

L; ■ bedded in the parent nebular regions, they are a unique labo- u i travio i e t (FUV), and near-ultraviolet (NUV) emission make a 

. _H , ratory for understanding the early dynamical evolution of the good combina t ion to study the properties of YSCs. For example, 

^ . associated stellar population, the properties of the Initial Mass ICorbelli et all d2009h have probed the stochasticity of the IMF 

»H ; Function (IMF), and the interaction between young massive stars using mu ltiwavelength photometry of a few tens of infrared se- 

, .0? , and the interstellar medium. Massive YSCs and complexes con- lected Y SCs in the inner disk of M33. In addition to the much 

tain hot massive stars which play an important role in galaxy re d UC ed distance ambiguities, the advantage of studying YSCs in 

evolution, by ionizing, heating and enriching the ISM of heavy external galaxies is that large samples can be easily selected to 

elements. They power the hot infrared (IR) radiation of dust and, probe different environments across the disk, such as spiral arms 

at the end of their lives, they eventually trigger further SF events and inte r-arms, central areas versus outskirts, and to compare 

dGrossi etalj2O10|>. The IMF is a probabilistic function describ- stal -.f orm ing sites in galaxies of different morphological types in 

ing the distribution in mass of stars at their birth; its logarith- a quiescent or dynamical stage 
mic slope above ~ 1 M has been often found close to the orig- 

inal value of -2.3 (Salpeter|[l955|). Given the rarity of massive The Local Group galaxy M33 is a quiescent blue galaxy at a 

stars predicted by the IMF, a statistically meaningful ensemble moderate inclination showing no sign of recent mergers or inter- 

of YSCs is needed to test the IMF at its upper end. action. Located at a distance D=840 kpc, it is currently forming 

External galaxies offer large numbers and a variety of SF stars at a rate of 0.45 M yr~'. The star formation rate (SFR) in- 

sites but, except for the Milky Way satellite galaxies, YSCs are ferred from the FUV emission is slightly higher (0.55 M yr _1 ) 

too compact for individual stars to be resolved. Even though the than that deduced by Ha emission (0.35 M Q yr~'), but the dif- 
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ference could be simply due to the uncertain correction for ex- 
tinction. When tr aced by the Ha emission (iKennicutd [l989t 
IVerlev et al.ll20 09). the SFR is found to decline radially with a 
scale length of about 2 kpc and then it drops sharply at galacto- 
centric distances of about 7 kpc. This location is often referred to 
as the star formation edge of the disk. This sharp drop is less ev- 
ident at other wavelengths, such as the mid-IR. In the far-UV 
the contribution of discrete sources drops at the SF edge but 
the diffuse emission has no apparent trunc ation or steeper de - 
cline out to the extent of the GALEX maps dThilker et al]l2005l) . 
The difference in the behavior traced by these indicators can be 
explained in terms of a SFR which has substantially decreased 
during the last 100-200 Myr, or by the lack of massive stars at 
large galactocentric radii. The former possibility has received 
some support by the recent d eep optical surveys (IGrossi et alj 
1201 it iDavidge & Puziall201 ll) of regions beyond the optical ra- 
dius where the average gas column density is only 2 x 10 20 cm -2 
dCorbellil 12003). These studies have revealed the presence of a 
population of relatively young main sequence stars with ages of 
about 100-200 Myr, indicating the existence of a more extended 
star forming disk in the recent past. Concentrations of ultraviolet 
light seen at large radii in some nearby spir al and dwarf galax- 



ies are often sign of i n situ star-formation dHunter et alj 201 1; 
iGil de Paz etaTT 120081: IThilker et all 120071) . Localized density 
enhancements beyong the bright star-forming disk may occur 
and give birth to a population of low-mass clusters. Clustering 
on sub-kiloparsec scale seen among stars with ages 10 Myrs 
close to the op tical edge of M33 has given evidence of recent 
star-formation (Davidge et al. 20lU). 

Inside the SF edge, the distribution of the ISM in the disk 
of M33 is not at all uniform. The flocculent spiral arms fade 
away around 4 kpc and this dramatic change might be re- 
sponsible for the variation in the properties of the interstellar 
clouds. Full imaging of the molecular cloud complexes in M33 
dEngargiola et all 120031: iHever etalj|2004t iGratier et all 1201 Oh 
has shown that those with masses a bove 1O 5 M are very rare 
beyond 4 kpc. iBastian et alj d2007l) also find that the number 
of young stellar groups also drops suddenly at radii > 4 kpc. 
However, a population of small mass molecular clouds exists at 
larger distances and it becomes dominant. This is also shown by 
the positive results of deep C O searches around d im IR sources 
close to the SF edge (see also Corbelr Tet al.ll2.01 Tb . 

Mid-IR maps are generally very helpful to select stellar clus- 
ters in their early evolutionary phase due to the presence of hot 
dust in the surrounding parent gas. The resol ution (6 arcm i n) and 
extent (16 kpc) of the 24 fim Spitzer map of IVerlev et al.1 (2007) 
is ideal for selecting YSC candidates throughout the whole disk 
of M33. The presence of UV or Ha emission is quite valuable to 
confirm the nature of mid-IR sources as YSC candidates since 
infrared emission alone cannot distinguish between YSCs and 
evolved stars, such as AGBs, that also contribute to the source 
counts in a galaxy (Cor belli et alj|201lb . The spatial extent of 
UV and Ha maps will confine our YSC analysis to somewhat 
smaller radii. 

Previous studies of infrared selected YSCs in M33 have been 
restricted to a few tens of them located in the inner disk. A more 
complete study of the properties of the whole YSC population 
throughout the star-forming disk and possibly beyond it, still 
need to be done. Finding YSCs in the outermost regions of M33 
will definitively prove the in situ origin of the pervasive stel- 
lar population of ages 100 Myrs recently discovered. The main 
aim of this paper is to investigate the properties of YSC popula- 
tion in M33 identified on the basis of their mid-IR, UV and Ha 
emission from the galaxy center out to about 12 kpc. By com- 



paring the emission of YSCs in the mid-IR, Ha, FUV and NUV 
bands with that predicted by theoretical models, we investigate 
possible differences in the distribution of mass, age, IMF and 
extinction throughout the M33 disk. This will help set further 
constraints on the spatial and time variation of the star formation 
process in this nearby blue galaxy. 

The paper is organized as follows. In Section 2 we present 
the sample and the multiwavelength dataset. In Section 3 we an- 
alyze the properties of the infrared and bolometric luminosity 
functions. The fitting method to the spectral energy distribution 
(SED) of YSCs is described in Section 4. In Section 5 we discuss 
the implications of the YSCs analysis on the IMF throughout the 
M33 disk. The cluster properties are summarized in Section 6 
and Section 7 lists our main conclusions. 

2. Multiwavelength data and the 24 fim catalog 

We select the young stellar cluster candidates of M33 using 
the 24/j.m images from Spitzer satellite data assembled by 
IVerlev etaD J2007). These authors cataloged 515 mid-infrared 
sources spanning over 4 orders of magnitude in luminosity and 
reaching, on the faint side, values equivalent to a single B2V 
star. However, this catalog is confined to the inner disk M33 
with only few sources extending as far as 7 kpc from the center. 
In order to extend the census of YSCs in the outermost regions 
of M33, we have us ed the largest avail able 24 fim map and the 
SExtractor software dBertin & Arnoutsl 19961) to identify d iscrete 
sources following the same method of IVerlev et all d2007l) . 

The sources were extracted above the 8<x level and their pho- 
tometry computed using the parameter "flux iso" in SExtractor, 
which uses isophotal photometry (sum of all pixels above thresh- 
old within the 8<x isophote); the background for the SExtractor 
photometry was the local median in a 32x32 pixel box centered 
on the source. We found 400 additional sources, so that our final 
list of 24 fim selected sources amounts to 9 1 5 objects (515 +400), 
whose spatial distribution is shown in Fig . [T| overlaid on the 
24 fim mosaic. It is evident that most of the IVerlev et alJ d2007l) 
sources (blue circles) lie in the inner region, whereas those iden- 
tified by us (red circles) occupy mainly the outer regions, extend- 
ing as far as 50 arcmin fr om the center . Table QJgives the position 
and size of the sources in lVerlev et al] d2007l) iden t ified a s YSCs. 
To the information already present in lVerlev et ail d2007l) we add 
the source size and its galactocentric radius. Table prefers to the 
newly added 400 sources; the reported size is the semi-major 
axis of the ellipse fitted to the 8<x isophote. 

In order to identify YSC candidates from the selected mid- 
IR sources, we complement the Spitzer 24 fim data with obser- 
vations in other bandpasses: 

i) Mid (8 fim) and far-infrared (70 and 160 /mi) data provided 
by the Spitzer InfraRed Camera (IRAC) and by th e Multiband 
Imagi ng Photometer for Spitzer, as described in IVerlev et"aT1 
d2007l) . 

ii) FUV and NUV maps o btained from the Galaxy Evolution 
Explorer (GALEX) satellite dGil de Paz et al.1 d2007l) . using the 
GALEX-UV images at their original sampling of 1.5 arc- 
sec/pixel. 

iii) The narrow-b and Ha map of Greenawalt et al. (1998), 
described in detail bv lHoopes & Walterbosl d2000h . 

We perform photometry in these other bandpasses using the 
position and size of the sources included in the 24 fim catalogue 
with the IRAF aperture photometry DAOPHOT package. Since 
the position of the source center and the extension of the emis- 
sion may change with wavelength, we have chosen a slightly 
larger photometric aperture, namely 1 .5 times the 24 fim size as 
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Fig. 1. The spatial distribution of 915 Mid-IR sources (circles) and 648 young stellar cluster s (squares) overlaid on the 24 fim 
Spitzer image. The light-blue and red circles in the on-line version represent IR sources from the lVerlev et al.l ([2007) catalogue and 
from the present study, respectively. 
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determined by SExtractor. We estimated the background (us- 
ing the "mode" function) in a 5-pixel wide circular annulus. 
Pixels 7cr below the background were rejected and the measure- 
me nts were corrected for foreground Galactic extinction follow- 
ing [SchlegiriLal] (DIM). 

2.1. The sample of Young Stellar Clusters 

The catalogue of YSCs of M33 includes all the Spitzer 24 yum 
sources with reliable photometry in the FUV, NUV, and Ha. We 
are limited by the smaller extent of the GALEX and Ha maps, 
compared to the mid-IR Spitzer maps, for the identification of 
24 /mi sources as YSCs. In the upper left panel of Fig.[2]we have 
marked the area covered by the GALEX maps together with 
the 24 /mi mosaic and all the 915 sample objects. The 24 fim 
map is larger than that covered by GALEX and 53 of the 24 yum 
sources lie outside of the latter. We were able to obtain reliable 
(error<0.5 mag) photometry for 789 sources in FUV, and 774 in 
NUV. A few other sources are too weak in the maps for an esti- 
mate _of_theJUy fluxes. The adopted zero-points in the AB sys- 
tem (lOkdl 19901) are 18.82 (FUV) and 20.08 (NUV). We added, 
in quadrature, to the formal error provided by DAOPHOT an un- 
certainty of 0.10 mag to allow for zero-point calibration errors 
(^ 0.05 mag) and large-scale background variance. 

The total field of view of the Ha map is 1.75 x 1.75 deg 2 
with a scale of 2 arcsec/pixel. The Ha map is also smaller 
than the 24 /mi mosaic (see Fig. [2] upper left panel) and 55 of 
the 24 /mi sources fall outside its borders. We obtained reli- 
able photometry (error<0.5 mag) for 772 sources. We added in 
quadrature an error of 5% in flux and corrected the Ha flux as- 
suming a 5% contaminatio n by Nil emission, as suggested by 
iHoopes & Walterbosl d2000h . 

IR sources with no, or weak, Ha counterpart are gener- 
ally faint and their nature is not obvious. A possibility is that 
they might be AGBs which are quite luminous in the mid IR. 
Alternatively, they could be small sites of star formation lacking 
massi ve stars, or deeply em bedded clusters, or background ob- 
jects. |McQuinnet^J l[2003) have assembled a catalog of AGBs 
in M33 based on Spitzer data. We have cross-checked our list 
with this catalog by searching for AGBs within the SExtractor 
24 /mi outer isophot. 

In conclusion, out of the original 915 24 /um sources, our fi- 
na l list of YSCs is m ade of 648 objects, 407 from the catalogue 
of lVerlevetal] d2007l) listed in Table [Q and 241 from the newly 
selected sources given in Table [2] The lack of the YSC label in 
the last column of Table [2] implies that either the MIR source 
is outside one of map boundary ( Ha, FUV or NUV) or it is a 
bona fide AGB (being associated to a MIR variable), or that it 
has unreliable photometry in Ha, or FUV- or NUV-band. We re- 
tain as YSCs the few sources associated to MIR selected AGBs 
which have a non-negligible number of ionizing photons, log 
Ha > 36.8 erg sec -1 : these are large HII regions with an AGB 
star from a previous generation-star or close by in projection. 

The YSCs selected at 24 fim have a wide range in size and 
luminosity and they are found in a variety of environments: from 
fully isolated, mainly in the galaxy outskirts, to crowded groups, 
especially in the inner spiral arms. Sources in Table|2]are mostly 
in the outer regions and isolated. In order to clearly assess the 
presence of UV and Ha counterparts and the reliability of the 
extraction process, we made up and inspected an atlas where, for 
each source, we assembled the images of the surrounding field in 
Ha, FUV, 8 and 24 /mi. This allowed us to select also a sample of 
24 /mi sources with insignificant FUV and Ha brightness as can- 
didate embedded star forming regions which will be examined 



in a forthcoming paper. An insight on their nature could come 
from a possible asso ciation with a molecular cloud or clump 
dCorbelli et al.ll201 lb but we have to wait for deep enough CO 
maps with sufficiently high spatial resolution to help in this re- 
gard. 

3. Radial densities and luminosity functions of IR 
selected sources 

Through out the paper we as sume a distance D to M33 of 
840 kpc dFreedman et al.llT99 lb to convert fluxes into luminosi- 
ties or to compute galactocentric distances. Assuming that the 
sources li e in the disk and adopting the geometry and warp 
model of ICorb elli & Schneider (1997), described in detail in 
ICorbelli & Saluccil d2000). we derive the distance Re from the 
center of our sample objects. The radial distribution of the 24 /mi 
sources and the selected YSCs are plotted in the top right panel 
of Fig. |2] We can see that the MIR sources are found up to the 
largest distances, while YSCs are distributed over the whole area 
of overlap of the various maps. 

In Fig [2] the bottom panels display the radial density distri- 
bution in linear and log scale, respectively. The densities of the 
24 /mi sources and of the YSCs behave similarly with an overall 
exponential decay with radius. In particular, we note i) a shallow 
slope in the inner disk, 0-4.5 kpc; ii) a steep decrease between 
4.5 and 8.5 kpc; and iii) a flattening in the outer disk. Past studies 
have found evidence for systemati c radial variat i ons in the star 
formation history of the M33 disk. lBastian et al.1 d2007l) showed 
that the number of young stellar groups drops suddenly at radii 
> 4 kpc and attributed this result to a change in the properties of 
the star forming sites. To study the inner and outer region clus- 
ters separately, we divide the 24 //m sources in 'inner clusters' 
with Re < 4 kpc and 'outer clusters' at larger distances. In this 
way, the size of the two groups is about the same. 

3.1. Estimate of the bolometric luminosity 

The bolometric luminosity of YSCs can be computed as the sum 
of the FUV and NUV luminosities uncorrected for absorption, 
the Ha luminosity corrected for extinction, L° Ha , and the total 
infrared luminosity L T ir. 

Lbo? ~ Lfuv + Lnuv + Ljir + 24-L° Ha (1) 

where the terms H a account for the ionizing radiation 
(lLeithereretal.il 19991) and Ljir for the UV radiation absorbed 
by grains and re-emitted in the IR. We have not considered the 
radiation longward of 2800 A which becomes important only for 
young clusters with luminosities lower than 10 38 erg s -1 . The lu- 
minosity in the FUV and NUV has been derived from 

Lw = vFyAnD 2 (2) 

with D the distance to M33, v = 1 .95 x 10 15 Hz for the FUV and 
v = 1.3 x 10 15 Hz for the NUV, and 

F v (Jy) = io 23 -(^+ 48 - 6 >/ 2 - 5 (3) 

AB indicates the corresponding magnitude in the FUV or NUV 
band. 

The total IR flux (F T ir) and L T ir are computed from the expres- 
sions: 



F TIR = 10' 14 x (19.5/^,(24) + 3.3F V (70) + 2.6/^(160)), (4) 
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Fig. 2. (Top left) Sky area covered by various maps along with the position of all the 24pm sources (915, dot symbols) and the 
sample of YSCs (648, open circles). Coordinates are for the J2000 epoch. (Top right) Histogram of the number of clusters N as a 
function of the galactocentric distance Rq for Mid-IR sources (white) and YSCs (grey). Bottom left & right panels Profiles of the 
number density of clusters per unit area in kpc 2 in linear and log scale, respectively, as a function of galactocentric distance for the 
same samples. See text for details on the disk deprojection model. 



Ltir = F nR x 10 3 xAnD 2 (5) 

where F TIR is in WnT 2 and F v in Jy dDale & Heloull2002h . The 
size of a large fraction of sources in our YSC sample is below 
the Spitzer telescope resolution limit at 70 and 160 pm. Hence, 
we cannot use the above expression to compute Ltir except for 
a limited sample of extended resolved sources. We can estimate 
Ltir from the emission at 8 and 24 pm if we find a good correla- 
tion between these and Ltir- For the sample of resolved sources 
at 70 and 160 pm we derive Ltir and its relationship to the 
emission at 8 and 24pm using the IRAC 8.0 pm photometry. 
This relation has been inferred previously but it might depend on 
the prope rties of the galaxy considered and on the scale exam - 
ined (e.g. lEngelbracht et al.l 120051: iPerez-Gonzalez et al.l l2006). 
To find a relation b etween the 8 and 24pm emission and Ltir, 
IVerlev et al.1 (H009) plotted the L24/LTIR ratio as a function of 
the 8/24 pm flux ratio for their sample of stellar complexes in 
the inner disk and derived a linear fit: 



log L Tm = log L 24 + 1 .08 + 0.5 1 log( -$-L) (6) 

where L24 is vF v x4nD 2 . We have checked that the above relation 
also holds in the outer disk i.e. that the numerical coefficients of 
the above relation for sources in the outer disk are, within the 
errors, compatible with those relative to sour ces in the inn e r disk . 
Thus, in this study we adopt the relation of [Verlev et al.l {2009) 
to compute Ltir from the 8 and 24pm fluxes. An error of 0.10 
mag is added in quadrature to the photometric error to compute 
the uncertainties in the 8 pm flux. The IRAC 8 pm map is smaller 
than the 24 pm mosaic (see Fig.[2]upper left panel) and 57 of the 
915 sources fall outside it. For these sources and for the few (25) 
without a reliable 8.0 pm photometry (error<0.5 mag) we use 
the relation 

logF v (8) = 1 .039 x log F v (24) + 0.002 (7) 

to derive the flux at 8pm which is needed to calculate Ltir- These 
coefficients are derived from a least square fit to the sources de- 
tected in both bands. 
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Extinction corrections to the Ha flux can be estimated using 
A Ha = 0.83Ay with A v 

A v = 1.7 xlog( L " R + 1) (8) 

LnUV + I^FUV 

This is sim ilar to the visual extinction formula used by 
ICalzettil J200I if the "fudge factor" used to infer the luminos- 
ity outside the FUV band is replaced with the luminosity in the 
NUV. 

3.2. Luminosity functions 

The luminosity function (LF) of YSCs is a fundamental probe 
of the star formation process of a galaxy. The cumulative LFs 
for the case of the inner and outer YSCs of M33 are shown in 
Fig- El We plot separately the results the cumulative function for 
the 24 fim flux for all the 915 24 yizm sources (Fig. a)), and 
the cumulative function for the TIR luminosity relative only to 
sources with an 8 /mi counterpart (Fig. [2b)). We see that the 
shape of the two LFs is similar; their slopes are given in Table 3. 
In the same figure we also plot the cumulative LFs of the selected 
648 YSCs as a function of Lt, i (Fig. [3jc)) and of the 24 /im flux 
(Fig. 12(f)). The values of the slope of the best least square fit to 
the data points between F v (24-fim) = 0.4 to 70 mJy are listed 
in Table 3 as 1 st fit, while the 2 nd fit includes points with higher 
values of F v (24fim). Since the values are lower than unity, the 
differential distribution N(L)dL oc L~ y dL has y lower than 2. 
Overall, the main features of the LFs are the following: 

(a) The slope of the cumulative function flattens at the faint 
end due to incompleteness. Crowding in the inner disk implies a 
somewhat higher completeness limit. 

(b) The LF of the inner disk shows a definite change of slope 
for a TIR or Bolometric luminosity of order of 10 40 erg s _1 . The 
brigh t end slope is simil ar to that found for brigh t HII regions in 
M33 dWvder et al.lll997t) and in late-type spirals (Caldwe lfet al.1 
11991b . 

(c) The LF of the inner and outer YSCs differ markedly: there 
are very few bright outer sources and the LF is steeper at the faint 
end. The presence of a break is not evident and, if present, it oc- 
curs at a lower luminosity. The distribution can be well described 
by a single power law. 

From the flattening of the distribution at the faint end, we 
estimate the completeness limit of the catalog around 0.4 mJy or 
Ljir ~ 5 x 10 37 erg s _1 . This co rresponds t o the bolometric lu- 
minosity of a sing le B1.5V star dCoxI feoOO). indicating that our 
sample is close to being complete even for faint obscured HII 
regions. The difference in the slopes of the inner and outer LFs 
between F v (24 yum) = 0.4 to 70 mJy can be explained by a dif- 
ference in the population of the star clusters. Namely, the inner 
regions have more massive clusters with the IMF fully popu- 
lated up to the upper mass end. Conversely, the outer regions 
form predominantly clusters of lower mass and hence the pres- 
ence of massive stars is rare and stochastic, as we shall discuss 
in Section 5. Owing to the similarity between the LFs for the in- 
frared and bolometric luminosity, we can exclude the possibility 
that YSCs in the outer regions are fainter in the IR due to the 
formation in an environment with a lower dust-to-gas mass ra- 
tio. As we will see below, the estimated YSC extinctions show 
no radial variation. In fact, beyond the edge of active star for- 
mation the MIR-to-FUV or to-Ha ratios increase. YSCs in the 
outer regions are intrinsically of lower luminosity: we will dis- 
cuss in the next Section whether this is due to an aging effect or 
to a different distribution in mass. 



As for the inner clusters, the distribution from the faint to 
the bright end can be broadly divided into two regimes with dif- 
ferent slopes. The steep slope at high luminosity is often ob- 
served in the LF of HII re gions, open clusters and associations 
dMcKee & Williams! 19971) . In the simple scenario, the change of 
slope represents the transition from poor to rich clusters, where 
the latter are populo us enough to reproduce the high-mass IMF 
(lOev & Cl arke 1998). In this case, the transition point between 
the two regimes marks the luminosity of the single brightest 
star: below this value, the observed statistics is modified by the 
sampling variance. We find that the transition occurs around 
F v (24fim) ~ 6 mJy, i.e. 10 40 er g s -1 , close to the luminosity 
of an 03V star (IVacca et alU 996). This implies that most of the 
bright MIR sources are in fact luminous YSCs. 

4. SED fitting 

Using the photometry in the FUV, NUV, Ha and MIR, we 
have constructed the SED of each Y SC. By a compari- 
son with model SEDs from Starburst99 dLeitherer et all 1 1999; 
IVazauez & Leitheredl2005l SB99), we can then derive the indi- 
vidual age, extinction, and mass; in addition, we provide some 
estimate of the average properties of the whole YSC population 
with implications on the IMF and on the ambient ISM. Below, 
we illustrate the procedure. 

4.1. Stellar population modeling 

We model the YSC SEDs with the SB99 single-age stellar pop- 
ulation synthesis. We consider that YSCs formed in an instan- 
taneous burst and as sume a sub-solar metallicity of Z = 0.004 
(Mag rini etall l2010h . The IMF is chosen to be the default 
Kroupa IMF, that is a broken power law (<p{M) oc M~ a ) with 
a slope a = 1.3 for . 1 M n < : M < 0.5 M and a = 2.3 for 
0.5 M < M < M max (Kroupal l200lh . Two sets of models were 
generated for different upper mass cutoffs: M max = 40 M and 
M max = 100 M , respectively. We use the Padova evolutionary 
tracks with full AGB evolution to account for their non negli- 
gible contribution to MIR fluxes. However, we did not consider 
the contribution of thermally pulsating AGB st ars that should 
be mi nimal for associations younger than 1 Gyr (Maraston et al. 
2006). The time evolution is followed with logarithmic sampling 
(30 steps) from 1 to 100 Myr. 

The SB99 calculations refer to a total stellar mass of 10 6 M Q , 
and no allowance is made for dust extinction or emission. The 
main output of the simulation is the emerging spectrum at each 
sampled age which we then convolve with the GALEX and 
Spitzer filter bandpasses to generate the synthetic photometry. 
The model spectrum we use is the one taking into account both 
stellar and nebular emission, the latter being computed for the 
case of ionization-bounded regions; under this hypothesis, SB99 
also provides the Ha luminosity. 

First, in order to include the effect of extinction and dust 
emission in the SB99 pho tometry, we select a dust mod el among 
3 cases made available by Weinga rtner & Draind (2001 PI: Milky 
Way with R v = 3.1 (MW31), LMC "average" distribution 
(LMC a vg), and LMC "2" distribution (LMC2). For our purposes, 
the difference among them is mainly in the relative extinction 
between FUV and NUV. We first assume that Ay is mainly pro- 
duced in the neutral ISM (screen geometry). Second, we choose 
a leakage factor as the fraction of the stellar radiation that 
escapes unimpeded by interaction with the ISM, that is without 

1 http://www.astro.princeton.edu/~draine/dust/dustmix.html 
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Fig. 3. The cumulative luminosity functions of: (a) the 24 /im sources (915), (b) the 8 fim sources (833) shown as a function of Ltir, 
(c) the selected YSCs (648) as a function of Lb„i and (d) the selected YSCs (648), as a function of the 24 jjm flux. Filled and open 
circles represent counts in the inner and outer regions, respectively. Linear fits for different flux ranges are shown as solid lines for 
the inner sample, and dashed for the outer one. 



any extinction and re-emission in the IR and also without ioniza- 
tion and Ha production. We consider two values of the leakage 
factor of and 1/3. Third, we allow the possibility that a fraction 
Kdust of the dust optical depth is contributed within the ionized 
phase (dusty HII regions): in this case, an important part of the 
ionizing flux is absorbed and reprocessed by dust, again without 
Ha emission. Again, we explore the influence of Kd us , for only 
two values, 0% and 30%, of the extinction within the HII region. 

In conclusion, given a dust model, a value of Ay, the leak- 
age fraction Fikg, and the fractional extinction in the HII region 
Kdust, we integrate the SB99 spectra and determine, apart from 
a multiplicative factor, the corrected FUV, NUV, Ha, and TIR 
luminosities, where TIR corresponds to the total amount of radi- 
ation absorbed by dust. 



4.2. x 2 minimization: Reddening, Age and Mass fitting 

As outlined above, the corrected FUV, NUV, Ha and TIR model 
luminosities depend on 3 cluster variables (age T, extinction Ay, 
and mass M) and 4 population variables (upper mass cutoff, dust 
model, leakage, and HII region dustiness). For a given set of 
values of the population variables, we used a^-minimization on 
each YSC and derived the best-fit values of the cluster variables 
(T,A V ,M) 



AT, A V , M) = 2 [lQ8(L obs )-loj(L model )-Af (Q) 

N obs 

where N is the number of bandpasses (4), L t, s is the observed 
luminosity, L mo d e i is the model luminosity, and cr 2 , is the esti- 
mated photometric error. In practice, the age and extinction of 
each cluster are determined by comparing the shape of the SED 
to that of the SB99 models, while the mass is the scale factor re- 
quired to best match the absolute values of the measured fluxes. 
Finally, A - log(M/ 1O 6 M ). As an illustration of the method, 
we display in Fig.|4]the SED best fits to four YSCs with differ- 
ent^ 2 values. 

The average {y 2 ) over the whole ensemble of YSCs of the 
minimized ^- 2 's is then used to discriminate among the different 
combinations of population variables. Such comparison is also 
performed separately for inner and outer clusters. Table |4] lists 
the mean and median (x 2 ) for all the combinations of cluster 
variables. For the inner clusters we find that (x 2 ) has a minimum 
for M max = 100, MW31, Fi kg = 0, and K dus1 = 0, while for outer 
clusters the minimum is obtained for M max = 100, LMC2, F^ g = 
0, and Kdusi = 0. The distribution of the^ 2 for such combinations 
is shown in Fig. [5] 

We conclude that the possibility of large-scale variations of 
the IMF, of radiation leakage, and of dusty HII regions do not ap- 
pear to be supported by these models that, however, favor a sys- 
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Fig. 4. Data (filled dots) and best fit model (open squares) SEDs for four YSCs with small errorbars but different^ 2 . Models in the 
bottom panels give good fits to the data. Models with ax 2 of order 100, as in the top-left panel, are of lower quality, but still yield 
acceptable fits. The few sources whose best fit model has^- 2 » 100 have a SED inconsistent with YSC models used in this paper. 



tematic variation in the dust model. The best fitting extinction 
curve for the inner and outer disk in fact differs, varying from 
Milky Way-type for the inner disk to LMC2-type for the outer 
disk. The best Of 2 ) is 14.8 for YSCs in the inner disk and 42.7 
for those in the outer disk. These values are rather large, despite 
the general goodness of the fit, because of our conservative un- 
certainty estimate. Models with (x 2 ) values lower than 17.1 and 
50.2 for the inner and outer disk, respectively, are within the 2<x 
probability level, if we normalize the^- 2 values to (x 2 ) = 14.8. 
Thus, they cannot be ruled out. This implies that some dust ab- 
sorption inside the HII region is still a possibility, but a 30% or 
higher fraction of leakage of radiation is unlikely. 

For each YSC, we can then compare the estimates of Ay and 
Lb i from the available photometry (see eqs. (1) and (8)) to the 
values obtained with the x 2 minimization method. This is illus- 
trated in Fig. [6] While the Lb„i values are in excellent agreement, 
the estimated A y& using eq. (8) are consistently higher than those 
inferred by the SED fits. The relation between the two quantities 
has a slope of 0.54. Since eq. (8) is purely empirical and does 
not take into account the evolutionary state of the cluster and the 
dust model, we consider the fitted value as more reliable. 

The age, mass and extinction distribution of the inner and 
outer YSCs as derived from the best-fit SED models are shown 
in Figure [7] Figure[8]and FigureHav respectively. 




Fig. 5. 

YSCs. 



log x 

Distribution of the x 2 for the sample of inner and outer 



5. The massive stellar population of YSCs 

Among all the models employed in the fitting, we find that the 
lowest average x 2 corresponds to a solution with the IMF pop- 
ulated up to 100 M , with negligible leakage of radiation and 
dust absorption inside the Stromgren sphere. From the (x 2 ) val- 
ues given in Table |4] we can see that an IMF with a maximum 



mass lower than 1OOM has a lower probability but it can- 
not be ruled out. Models with an IMF fully populated only up 
to 40 M have been used also to fit the SED of clusters with 
Lboi < 2 x 10 39 ergs~' only, or the sample of clusters beyond 
7 kpc (the nominal edge of the star-forming disk). However, 
this choice did not provide an improvement in the average (x 2 )- 
Therefore, we consider models with the upper mass cutoff at 
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Fig. 6. (Left panel) Comparison of the Ay values for the YSC sample derived from two different methods: Ay as given in eq.(8) 
using the TIR luminosity (A V ,R ), and using the SED model fits (Ay'). (Right panel) Same as the left panel, but for the luminosity 

estimated from the sum at different wavelengths an d from the SED model fits In both cases, the best linear fit is shown 

by the solid line. 
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Fig. 7. Histogram showing the distribution of YSC ages in the inner and outer disk of M33. On the x- axis the age 'T is in a log 
scale; on the y- axis the number of clusters for each age bin is normalized to the total number of clusters in the inner (left panel) and 
outer (right panel) disk. 



1OOM as the best representation of the stellar population of 
YSCs throughout the M33 disk. 

There is a clear difference between the mass distribution 
of clusters which form in the inner disk with respect to those 
forming at large radii: beyond 4 kpc clusters more massive than 
1O 4 M are very rare. Similarly, beyond the Ha edge at 7 kpc, 
most YSCs have masses smaller than 1000 M Q . If clusters form 
with total mass below 1O 4 M , the number of stars is not high 
enough for the IMF to be fully populated up to 100 M Q . It 
is well known that the most massive stars in a cluste r tend 
to have a mass that increases with the cluster mass dLarsonl 



Il982l) . This correlation leads to the question of whether the up- 
per limit to the stellar mass in a particular cluster depends on 
the cluster mass because of some physically limiting process, 
or it is simply the result of random sampling of the IMF. In 
the first case, low-mass clusters cannot produce high-mass stars 
dWeidner & Kr oupa 2004, 2006). In the second case they could 
as long as there is enough gas, and intermediate-mass clusters 
should occasionally be found with unusually massive stars - 
"outli ers" in the IMF (Elmegreen 200(1 ICervino & Luridianal 
2006). A consequence of this model, often referred to as stochas- 
tic or randomly sampled IMF model, is that the summed IMF 
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Fig. 8. Distribution of the cluster mass in the inner and outer disk of M33. We plot the log of mass (in solar units) and of N+l 
(where N is the number of clusters in each mass bin), respectively. The grey and white histograms are for the masses computed with 
or without incompletness corrections (see Section 5), respectively. 
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Fig. 9. Distribution of the visual extinction in the inner and outer disk of M33. Ay is in magnitudes, while the number of clusters 
per bin is shown normalized to the total number in each sample. 



from many small mass YSCs should be the same as the IMF of 
a massive YSC. Considering the case for M33, the higher (x 2 ) 
shown by the IMF model with a mass cutoff lower than 100 M , 
even considering only the low-luminosity cluster sample or clus- 
ters beyond a certain radius, does not favor models which relate 
the maximum stellar mass to the cluster mass. We recall that 
the res ults on the clus ter birthline of M33 support the stochastic 
model (Corbelli et al. 2009), which also app ears favoured by the 
analysis of the IMF in more distant galaxies dCalzetti et alj2 010; 
Fumagalli et al. 201 1). Hence, in what follows we will compare 
our data to the predictions of the stochastic IMF. 

We simulate the cluster light and mass distributions consid- 
ering a stochastic IMF from 0.1 to 100 M with slope a = -2.3 



for stellar masses greater than 0.1 M and a = -1.3 for lower 
masses. We simulate 40,000 clusters that are distributed in num- 
ber according to their mass between 20 and 10 4 M Q as 

N(M)dM = KM s dM (10) 

where K is a constant, and 5 is the spectral index of the Initial 
Cluster Mass Function, which we assume to be 5 - -2. This 
is in close agreement with the mass distribution of our IR se- 
lecte d YSCs (see also next Sectio n) and with previous findings 
(e.g. Ide Griis & Parm entier 2007, and references therein), even 
though the results of the simulations are not very sensitive to the 
value of 6. We populate each cluster with stars randomly selected 
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from the st ellar IMF. More detai ls on the numerical simulation 
are given in Corbel li et a i] (l2009h . 

Assuming the mean value (x 2 ) = 14.8 of the best fit model 
for the inner disk YSC distribution as the peak of the parent dis- 
tribution, we consider YSCs whose best fit model has (^ 2 ) < 100 
well fitted by the SED models. This is shown also by Figure [4] 
Leaving out clusters with^ 2 > 100, with large photometric er- 
rors (greater than 2 in log) and clusters older than 25 Myr, we 
now compare the relation between the FUV and Ha luminosi- 
ties with that obtained by simulating a distribution of clusters at 
birth. We correct the observed FUV and Ha luminosities for ex- 
tinction using the extinction value of the best fitting SED model. 

In the bottom panel of Fig. [10]we show the distribution of Ho- 
luminosities as a function of the FUV luminosities corrected for 
extinction according to the best fitting SED model. The dashed 
line shows the linear relation which is consistent with the dis- 
tribution of the bright clusters. Clearly, the distribution deviates 
from the linear trend as Lfuv drops below 10 39 erg s _1 . The sen- 
sitivity of our survey is good enough to trace the non-linear rela- 
tion for Lfuv > 5 x 10 37 erg s _1 . In the upper panel of the figure 
the same data are shown together with the average value of log 
Liia (filled squares) in 0.4 wide bins of logLfw The filled cir- 
cles show the same expected values as derived from the simula- 
tion of the stochastic IMF for clusters at birth (i.e. younger than 
3 Myr). The most relevant result is that the expected deviations 
from the linear scaling for a random sampled IMF model agree 
with the data. The measured Ha luminosities for bright clusters 
are about 0.2 dex lower than those predicted by the stochastic 
model. This is because our sample does not contain very young 
massive stellar clusters, and aging for massive clusters lowers 
the expected Ha luminosity even if T< 10 Myrs. From the fig- 
ure we also see that stochasticity implies a large spread in the 
Ha luminosities for L^uv < 10 39 erg s _1 . 

The presence of outliers in low-mass clusters populated ac- 
cording to a stochastic IMF implies that the relation between the 
cluster luminosity and its mass is no longer linear, or uniquely 
determined. We show in Fig.QT]the distribution of cluster masses 
as a function of the cluster bolometric luminosity at birth. We 
have simulated the mass distribution using for massive clusters a 
M - Li,oi relation consistent with that given by the SED models 
for objects younger than 3 Myr. The filled dots in the Fig. [TT] are 
the average values of logM, in solar units, for different values of 
log Lb„i. For low-mass clusters, stochastic sampling of the IMF 
implies a mass higher, on average, than that obtained from the 
extrapolation of a fully populated IMF. Deviations from a linear 
scaling are evident as L\, \ approaches 10 40 ergs _1 and become 
relevant for the average mass when L bo i < 10 39 erg s _1 . A clus- 
ter of given luminosity can be of smaller mass when some bright 
outlier dominates its integrated light, or can have a larger mass if 
no outliers are present. The dashed line in the Fig. QT] shows the 
linear approximation to the logM - logL\, i relation for clusters 
at birth. We use this relation to derive the cluster mass from the 
fitted bolometric luminosity. Given a cluster age T, we define a 
critical luminosity 

logL%l = 413-l.2logT (11) 

and, when Lj, i < L c h "' r we correct for incompleteness using the 
expression: 

log^- = 0.58logL bo i - 24.45 + O.llogT. (12) 
M 

In the above formulae Lboi is in erg s _1 , T is the cluster age in yr 
for clusters older than 3 Myr and is T — 10 6 yr for younger clus- 
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Fig. 10. (Bottom panel) Distribution of Ho- luminosities as a 
function of FUV luminosities corrected for extinction accord- 
ing to the best fitting SED model. Units are erg s _1 . The dashed 
line shows the linear relation consistent with the distribution for 
bright clusters. Filled triangles and open circles are for inner 
and outer disk YSCs, respectively. (Top panel) The same data 
are shown using small dots together with the average value of 
logLua (filled squares) in 0.4 wide bins of logLfuv- The filled 
circles show the expected values as derived from the simulation 
of the stochastic IMF for clusters at birth. 

ters for which the time evolution of L\, i is negligible. For Lb i > 
U. n \ the SED models show that the cluster mass is well approxi- 
mated by the relation logM/M^ - logL bo i - 44.3 + 1 2logT . 

Our main conclusion of the analysis of YSCs in M33 is that 
a stochastically sampled IMF with a maximum mass of 1OOM 
well describes many properties of their integrated light distribu- 
tion throughout the whole disk. 

6. YSC properties 

We have corrected the mass of the less luminous clusters inferred 
by the SED model fits for incompleteness of the IMF according 
to the results of the previous Section. The average mass of low 
luminosity clusters is underestimated if no correction is applied. 

The mass-radius relation of the YSCs of our sample is illus- 
trated in the right panel of Fig. [12] Here, the size is the 24jum 
semi-major axis from SExtractor. The distribution shows a large 
sc atter for sizes less than 10 pc, as alre ady noted for our Galaxy 
bv lZepf et al.l (1991) and lLarsenl (|2004) who found a poor mass- 
radius correlation for compact star clusters (radius < lOpc) with 
ages between a few tens to hundreds of Myr. In nearby spiral 
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Fig. 11. Simulated distribution of cluster mass as a function of 
the cluster bolometric luminosities at birth. The large filled cir- 
cles show the average cluster mass in each bolometric luminosity 
bin. The broken dashed line is a fit to the average logM - logL\,„\ 
relation: at high luminosities it follows the linear relation be- 
tween mass and luminosity expected for a fully populated IMF, 
while for log L/, / < 3 x 10 39 erg s _1 it approximates well the 
simulation results. 



galaxies, iLarsenl £2004) obtains M oc R 21 in relatively young 
comp act clusters wi th age < 1 Gyr and R < lOpc (confirmed 
bv IGrossi etai1l2010l for YSC in M33). In our sample YSCs 
larger than 10 pc show a good correlation with a log-log slope 
of 2.0 9 + 0.01 (least square fit). In M33 iBlitz & Rosolowskvl 
(l2005h obtain a similar slope considering giant molecular clouds 
(M c i oc R 22 ), simi lar to other nearby galaxies such as M51 
dBastian et al.ll2005l) for GMCs with similar ages than the YSCs 
in M33 but larger sizes. We may thus arrive at the conclusion 
that the YSCs display average characteristics reminiscent of the 
protocluster environment from w hich they we r e born . A simi- 
lar conclusion was also drawn by IGrossi et alJ (1201 Oh from the 
analysis of stellar clusters and complexes. 

The relation between mass and Ha luminosity is shown in 
the right panel of Fig. [12] We observe a good correlation for 
cluster masses in excess of ~ 500 M Q , although the spread is 
significant at the faint end, as expected for the stochastic IMF 
model (see previous Section). 

The age distribution of the inner and outer YSCs is dis- 
played in Fig. [7] The distribution of the former is narrower 
than the latter and indicate that clusters tend to be younger 
in the inner regions of M33. On the other hand, the cluster 
mass distribution for the two subsamples shown in Fig. [8] re- 
veals a broader spread for the inner sample, while the outer 
regions peak at ~ 1O 2 M and are devoid of massive objects 
(> 10 4 M Q ). The distribution of cluster masses, corrected for in- 
completeness of the IMF, is compatible with a differential dis- 
tribution of index -2. This is consistent with the value found 
in the MW and in M33 Jyan den Bergh & Lafontaind 1 19841: 
IChristian & Schommerlll988Llde Griis & Parmentieril2007l) . 



For the large majority of the YSCs selected in the mid-IR 
with a UV and Ha counterpart we find ages between 3 and 
lOMyr. Given the sample selection criteria, we are unable to 
trace the very early phases of the cluster formation process. 
Beyond the Ha edge, at about 7kpc, we do not find clusters 
younger than 6Myr and the average age is of order lOMyr (see 
Figure fT3l). There is no clear evidence that this result is due to 
a selection bias as if, for example, outer YSCs, being of smaller 
mass (and luminosity), would be more difficult to detect in their 
early stages, especially due to the requirement of an Ha and UV 
counterpart. The lack of younger clusters in the outer regions is 
not ameliorated if the IMF is truncated at 40 M , while this IMF 
model increases slightly the number of clusters with T< 3 Myrs 
in the inner regions. 

Fig. [9] shows the Ay distributions. Most of the YSCs in our 
selected sample have Ay between 0.2 and 1. We do not find ev- 
idence that the outer clusters suffer from lower extinction, al- 
though there is evidence for a different extinction curve. The ra- 
dial variation of age and cluster mass for our sample of YSCs is 
shown in Fig. [13] We can clearly see that, with increasing galac- 
tocentric radii, the clusters become less massive and older. The 
SEDs of clusters at very large radii are characterized by a higher 
IR-to-UV and IR-to-Ha light ratios. This property is well fitted 
by the synthesis models of older clusters. However, the pres- 
ence of low-mass clusters as young as ~10Myr beyond the SF 
edge of the disk, where the average HI column disk is only a few 
times 10 20 crrT 2 , is remarkable. Localized density enhancements 
are evidently present with the consequent formation of small H2 
clouds, birthplaces of low mass YSCs. The molecular gas disper- 
sal tim e is comparable to the YSC ages that we find (Leisawitz 
1990). Hence only some residual gas is present in the YSC en- 
vironment, similar to that detected around MIR sources wit h 
Ha counterparts close to the SF disk edge dCorbelli et alj|201lb . 
Hence, the low extinction values found are due to our selection 
criteria: we have included in our sample only MIR sources with 
FUV and Ha counterparts. 

Finally, in Fig. [14] the relations between age, mass and Ay 
are displayed as density images to overcome the crowding of 
points. Again, we can see that the distributions are rather tight 
with most YSCs having Ay =0-1 mag, log(T) =6.5-7.1 Myr and 
log{M/M Q ) =2-4. The peak of the density maps corresponds to 
A v = 0.3 mag, log(T) = 6.7 Myr, and log(M/M ) = 2.5. 

7. Summary and Conclusions 

In this study we have selected 915 mid-IR sources in M33 from 
the largest available 24 /mi Spitzer map. Complementing the mid 
and far-infrared Spitzer data with UV data from the GALEX 
satellite and with Ha data, we have investigated the properties of 
648 mid-IR sources which are YSC candidates distributed from 
small to large galactocentric radii. The observed SED for each 
YSC is compared with theoretical synthesis models to derive 
the cluster age, mass and extinction. SED models with differ- 
ent leakage fractions, IMF upper mass cutoffs, extinction curves 
and dust absorption fraction in HII regions have been tested sep- 
arately for the inner (Re < 4 kpc) and outer cluster sample. This 
allows us to constrain the IMF and ISM properties of the inner 
and outer disk of M33. The main results can be summarized as 
follows. 

- We find IR sources as far as 16 kpc from the center, corre- 
sponding to the boundary of the 24 fim mosaic and to the 
extent of the warped HI disk. The completeness of the MIR 
source catalogue presented here is about 0.4 mJy, equivalent 
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Fig. 12. (Left panel) Distribution of cluster mass as a function of size. The best fit for clusters with size > lOpc is shown by a solid 
line. Masses have been corrected for incompleteness. (Right panel) Relation between cluster mass and the Ha luminosity. 
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Fig. 13. The radial distribution of mass, age and L TIR /L FU v of YSCs as a function of the galactocentric distance R G - Circles indicate 
YSCs with unreliable fits Of 2 > 100). 
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Fig. 14. Iso-density contours showing the age-mass-reddening relationship for the YSCs in our sample. 
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to the luminosity of a single B1.5V star. The radial decline 
of the IR source density becomes steeper at galactocentric 
distances R c =4-5 kpc and flattens out beyond the optical ra- 
dius (8.5 kpc). YSCs selected in the mid-IR with Ha and UV 
counterpart follow a similar radial distribution. 

- The luminosity distribution of 24 fim sources as a function of 
the infrared or bolometric luminosity is steeper at the faint 
end in the outer regions, which are devoid of very bright 
sources. The LF in the inner disk displays a double slope, 
steepening at the high luminosity end where the slope is sim- 
ilar to that observed for HII regions in M33. 

- The average extinction of the YSCs is modest, since most 
of them have Ay < 1. The similarity between the bolomet- 
ric and infrared LFs throughout the M33 disk and the lack 
of radial trends in A v imply that the dusty environment does 
not vary appreciably between star forming sites at small and 
large galactocentric radii. However, a Milky Way-type ex- 
tinction curve is preferred in the inner disk while for the outer 
disk an LMC2-type gives the best fitting SED models. 

- The mass-radius relation of YSCs shows a linear correlation 
with slope 2.09 + 0.01, in close agreement with the GMC 
mass-size relation. Hence, the stellar densities of the YSCs 
in our sample are representative of the protocluster environ- 
ment. 

- Low luminosity clusters (M < 1O 3 M ) become the domi- 
nant population at large galactocentric radii. We find YSCs 
with age of about lOMyr as far out as 12 kpc, correspond- 
ing to the extent of the Ha map used for the identification of 
YSC candidates. Previous surveys have uncovered a popula- 
tion of ~100Myr old stars in the outskirts of M33, but our 
finding implies a more recent episode of star formation in the 
outer galaxy. 

- The YSCs in the inner disk show a broad mass and age range, 
extending in mass up to 10 5 M G , and with ages between 3 
and lOMyr. The peak of the YSC distribution lies at Ay = 
0.3 mag, age=5 Myr and mass= 300 M . 

- The best fit synthesis models to the SED favour an IMF with 
an invariant upper mass cutoff at 100 M Q throughout the M33 
disk. Taking into account the stochasticity of the IMF, rele- 
vant for low-mass clusters which cannot populate the IMF 
up to the upper mass end, we are able to reproduce the devi- 
ations from linearity of the UV-continuum-Ha relation that 
we observe for YSCs. 
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Table 1. The position, size and the galactocentric distances 'Rc' for sources in IVerlev et~a"fl (120071) classified as Young Stellar 
Clusters and included in the analysis of YSC properties presented in this paper. The complete table is available in the electronic 
form only. 



ID 


«(2000) 


5(2000) 


Size 


Rg 




(degree) 


(degree) 


(arcsec) 


(kpc) 


1 


23.374977 


30.309910 


15.3 


5.37 


2 


23.361427 


30.290756 


3.8 


5.66 


3 


23.372110 


30.301689 


7.1 


5.49 


4 


23.332920 


30.305302 


5.1 


5.47 


5 


23.335817 


30.304033 


3.6 


5.48 


514 


23.599760 


30.950455 


7.5 


4.60 



Table 2. Additional 24 fim source catalogue (added to the catalogue o flVerlev et al.l(l2007l) ). The position, size and the galactocentric 
distances '/?c' are given in addition to fluxes at 8 and 24 /mi. The '-' symbol in column 7 and 8 indicate sources outside the 8 /urn map 
boundary. In the last column the symbol "YSC" indicate sources classified as Young Stellar Clusters and included in the analysis of 
YSC properties presented in this paper. The complete table is available in the electronic form only. 



ID 


«(2000) 


5(2000) 


Size 


F 24 yum 




F 8 yum 


cr 


Rg 


Note 




(degree) 


(degree) 


(arcsec) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(kpc) 




516 


23.320433 


29.687447 


1.8 


1.072 


0.020 






15.31 




517 


23.397005 


29.691725 


1.6 


0.850 


0.016 






14.85 




518 


23.264629 


29.772163 


3.8 


15.433 


0.091 






14.53 




519 


23.341816 


29.774063 


1.2 


0.310 


0.014 






14.02 




520 


23.529587 


29.812027 


1.9 


1.118 


0.016 






12.60 




521 


23.348316 


29.812914 


1.8 


1.011 


0.017 






13.45 




783 


23.457613 


30.742470 


3.5 


1.779 


0.076 


3.918 


0.223 


1.35 


YSC 


784 


23.308142 


30.744352 


3.5 


1.280 


0.057 


1.384 


0.109 


3.69 


YSC 


785 


23.450867 


30.746641 


8.0 


17.984 


0.317 


36.865 


0.676 


1.47 


YSC 


786 


23.823700 


30.746367 


1.1 


0.277 


0.007 


0.096 


0.014 


7.78 


YSC 


915 


23.571716 


31.531706 


1.4 


0.336 


0.013 






13.74 





Table 3. The slope of best least-square fit to the cumulative distribution of the clusters of inner and outer region of M33 galaxy. The 
sample 'a' (24yum sources (915)) and 'd' (selected YSCs (648)) represents the cumulative distribution of the clusters as a function 
of 24 /im flux whereas the sample 'b' (8 //m sources (833)) and 'c' (selected YSCs (648)) represents the cumulative distribution 
of the clusters as a function of Ljm and Lb„i respectively. The 1 st and 2 nd fit represents the fitting in the lower and higher flux (or 
luminosity) ranges. 



Sample 


Region 


1" fit 


2"" fit 






slope ±<t 


slope ±<x 


(a) 


inner 


-0.44 ± 0.02 


-1.12 + 0.06 




outer 


-0.84 ± 0.03 




(b) 


inner 


-0.42 ± 0.02 


-1.38 ± 0.11 




outer 


-0.88 ± 0.03 




(c) 


inner 


-0.32 ± 0.02 


-1.40 ±0.07 




outer 


-0.68 ± 0.04 




(d) 


inner 


-0.38 + 0.02 


-1.12 ±0.06 




outer 


-0.78 ± 0.02 
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Table 4. The sum of the^ 2 obtained for different set of model parameters for inner and outer disk clusters. 



IMF mass- 


EXT curve (MW/ 


K-dust 


Fikg 


X 1 Inner 


X 1 Outer 


cut off (Mq) 


LMCW/LMC2) 


(%) 


(%)) 


(Mean) 


(Mean) 


100 


MW 








14.8 


44.5 


100 


IMC AV 








18.6 


45.9 


100 


LMC2 








17.3 


42.6 


100 


MW 


30 





15.7 


49.2 


100 


IMC AV 


30 





21.3 


53.5 


100 


LMC2 


30 





18.6 


49.0 


100 


MW 





30 


21.1 


68.4 


100 


LMCav 





30 


21.3 


69.2 


100 


LMC2 





30 


20.3 


67.9 


100 


MW 


30 


30 


22.3 


66.2 


100 


LMC AV 


30 


30 


24.4 


70.3 


100 


LMC2 


30 


30 


21.4 


66.8 


40 


MW 








15.4 


50.7 


40 


IMC AV 








19.0 


52.0 


40 


LMC2 








17.3 


48.0 


40 


MW 


30 





19.5 


58.9 


40 


LMC AV / 


30 





25.8 


65.2 


40 


LMC2 


30 





21.0 


58.6 


40 


MW 





30 


29.1 


82.4 


40 


LMCav 





30 


27.3 


82.2 


40 


LMC2 





30 


26.4 


80.9 


40 


MW 


30 


30 


46.9 


94.6 


40 


LMCav 


30 


30 


52.5 


102.1 


40 


LMC2 


30 


30 


45.6 


95.7 
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